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We present for the first time to the nuclear physics community the Hamiltonian Mean 
^ ■ Field (HMF) model. The model can be solved analytically in the canonical ensemble 
OO . and shows a second-order phase transition in the thermodynamic limit. Numerical mi- 
Q ! crocanonical simulations show interesting features in the out-of-equilibrium regime: in 
t^^ I particular the model has a negative specific heat. The potential relevance for nuclear 

multifragmentation is discussed. 
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^ , 1. INTRODUCTION 

C ■ 

In the last years there has been a lot of interest in the nuclear physics community 

for multifragmentation reactions and for the study of liquid-gas phase transition in finite 

rS \ systems such as nuclei [0-|^ . It is still a debated question whether nuclear multifragmenta- 

c^ ■ tion is an equilibrium type of transition, and if it is first or second-order |]7|-|15| . With the 



scope of adding new useful arguments to this debate, we present for the first time to the 
nuclear physics community the Hamiltonian Mean Field (HMF) model, a system of classi- 
cal spins coupled through long-range attractive forces. HMF can be solved analitically in 
the canonical ensemble and shows a second-order phase transition in the thermodynamic 
limit. Moreover HMF has the double advantage of allowing a microcanonical and dynam- 
ical approach to explore the dynamics of a phase transition in a finite system. We believe 
that, in the same spirit of the Ising model, introduced many years ago, but still stud- 
ied with extreme interest in statistical mechanics, the HMF model, though its dynamics 
is probably less complex than that one of a real system, can become a tool of primary 
importance to extract information on the microscopic dynamics of a second-order phase 
transition. 

In this paper we discuss the equilibrium and out-of equilibrium properties of the model, 
focusing our attention on the importance of the relaxation to equilibrium and on the 
possible discrepancies between different ensembles. In this respect the model can be 
extremely useful also for the nuclear multifragmentation phase transition. In particular 
in section 3, we show how misleading information on the order of a phase transition can 
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Figure 1. Equilibrium microcanonical simulations for N=500 (open squares) in compari- 
son with the theoretical predictions (full curves). The numerical results are averaged over 
several runs. The small discrepancies around the critical point {Uc = 0.75) and in the 
high energy region are due to finite-size effects. 



be obtained from a specific heat analysis [|13| if the system is not perfectly equilibrated. 



2. HMF: THE CANONICAL ENSEMBLE 

The Hamiltonian Mean Field (HMF) model was introduced by Ruffo [0 in 1994 and 
has been since then intensively studied both analitically and numerically [p!7|-|21|. The 
Hamiltonian of the model describes a system of N fully-coupled particles moving on the 
unitary circle, where the particles are characterized by the angles di and the conjugate 
momenta pi. It reads: 
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being K and V the kinetic and potential energy. If one considers a spin vector associated 
to each particle rrij = [cos6i, sinOi], the Hamiltonian then describes a linear chain of 
A^ classical fully-coupled spins, similarly to the XY model, and we can define a total 
magnetization vector M = [Mx, My] = -^ J2iLi mj • A very important feature of the HMF 
model is that it presents an analytical solution in the canonical ensemble. The solution has 
been worked out by Antoni and Ruffo, by means of the Hubbard-Stratonovich method 



to perform the integration over the angles ||T6|. The system behaves as a ferromagnet 



at low energy and shows a second-order phase transition at the critical energy density 
Uc = 0.75 {U = E/N where E is the total energy), corresponding to a critical temperature 
Tc = 0.5 p6|-p!8|. If the sign of the interaction is inverted, one gets an antiferromagnetic 



behavior: this case has been studied in detail in refs []I6|,^^. The theoretical caloric 
curve reads: 

N dp (2/3) ^ 2 V ^ ' ^ ^ 

where (3 = 1/kT is the inverse of the temperature (the Boltzmann constant k is set 
to be equal to 1). The canonical caloric curve is shown in fig.l (a) as a full curve: it 
increases almost continuously with energy, but has a discontinuity at the critical point. 
The modulus of the total magnetization, M, shown as full curve in fig. 1(b), is close to one 
for very low energy and decreases to zero by increasing the energy. It is equal to zero for 
U > Uc- It has been shown that the phase-transition in the HMF model is a second-order 



one with mean field critical exponents f3 = 1/2, a = Q |16,18]. It is very intriguing the 
fact that, though HMF was not constructed to describe a nucleus, its caloric curve turns 
out to be very similar to the one measured by the ALADIN group [|l| and more recently 
by other experimental collaborations |^,^,^ for nuclear systems. 

3. HMF: THE MICROCANONICAL ENSEMBLE 

In this section we discuss the microcanonical numerical simulations and we compare 
with the equilibrium canonical solution. The Hamilton equations of motion for the A^ 
particles are 

0^=p^ , Pi = -sin{9i)M, + cos{9i)My , i = l,...,N , (3) 

and can be solved numerically for different sizes of the system. In this way HMF allows a 
dynamical microcanonical approach, perfect to explore the dynamics of a phase transition 
in a finite system. The equations of motion were integrated by means of a 4th order 



simplectic algorithm ||2^ with a time step 6t = 0.2 and a relative error in the total 
conserved energy smaller than AE/E = 10~^. More details can be found in Refs |17,18 



3.1. Equilibrium Results 

In this paper we consider a system with A^ = 500, a number of particles of the same 
order of the typical number of nucleons in a nuclear multifragmentation reaction. We start 
the system out-of-equilibrium, in the so-called "water bag", i.e. by putting all rotators 
at gj = and giving an initial velocity according to a constant probability distribution 



function of finite width centered around zero |16,17]. We remind the reader that, though 
the interaction in HMF is not explicitly constructed to model a nuclear system, such a 
far-off-equilibrium initial condition is a good way to simulate the strong off-equilibrium 
conditions present in a hot and compressed nuclear system before multifragmentation. 

We follow the dynamics of HMF up to a final time t = 2 ■ 10^ and we show in fig.l the 
final results of the numerical integration. The numerical results were integrated also on 
several different runs. Such long integration times are necessary because, when out-of- 
equilibrium initial conditions are used, the system equilibration can be very slow. This is 
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Figure 2. We show, for a typical event, the time evolution of temperature T for A^ = 500 
and different energy densities U = 0.4, 0.6, 0.7, 0.8. The dashed lines are the corresponding 
canonical equilibrium values. Note that for U=0.8 the dynamical temperature converges 
to a slightly higher value (dotted line), which is the equilibrium one considering finite-size 
effects, see text. 



especially evident close to the critical point, where the presence of quasi-stationary (long 
living) non-equilibrium states (QSS)[|r^,|T^ has been revealed. The temperature, the total 
magnetization, the kinetic fluctuations and the specific heat per particle are reported as a 
function of the energy density as open squares and compared to the theoretical predictions 
(full curves). The temperature is computed from the average kinetic energy per particle 
T = 2 < K > /N, where the symbol < > stands for time averages. The kinetic energy 
fluctuations S, obtained from the scaled variance of the kinetic energy E = ^^^ /^^^^ > 
are compared to the theoretical prediction in the microcanonical ensemble. The latter is 

given by the formula E = ^y 1 — [1 — 2M(^)] [|^. Finally, in panel (d), we compute 
the specific heat per particle at constant volume in the microcanonical ensemble using 



the formula proposed for the first time by Lebowitz, Percus and Verlet [^ 
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In general, the numerical microcanonical caloric curve and the magnetization obtained 
at very long times agree quite well with the canonical predictions, apart from small finite- 
size effects located mainly around the critical point and in the supercritical region, the 
homogeneous phase. In the latter, the magnetization, due to the fact that N=500, is not 
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Figure 3. Non-equilibrium microcanonical simulations N=500 (open circles) in compar- 
ison with the theoretical predictions (full curves). The integration was stopped at time 
t = 3000 and the results averaged over several runs. Compare also with figs. 2,4 and see 
text for further details. 



exactly zero and assumes a finite value, i.e. M ~ 1/yN. This has visible consequences on 
the temperature, according to formula (2) and thus also on the kinetic energy fiuctuations 
and the specific heat. We note that the microcanonical specific heat at equilibrium is 
always positive and agrees with the theoretical formula (reported as a solid line). In 
conclusions the canonical ensemble and the microcanonical simulations at equilibrium are 
in good agreement. The equivalence of different ensembles in the thermodynamic limit is 
also supported by recent analytical results in the microcanonical ensemble |]2l| . 



3.2. Out-of-Equilibrium results 

In the previous subsection we have presented the microcanonical numerical results ob- 
tained in the equilibrium regime. In the present one we study the process of relaxation 
to equilibrium and we show that HMF has a very rich dynamics. 

In fig. 2 we report the time evolution of the temperature for different energies, and we 
study how it relaxes to the canonical value, represented by dashed straight lines. Though 
for U=0.8,0.4 the equilibrium temperature is almost reached already at a time t = 3000, 
close to the critical point, and in particular for the energy U=0.7 , we have to wait for 
much longer times (for a time t > 5 • 10^). In fact, before the relaxation to the canonical 
value, there is a well defined plateau, where the temperature assumes a constant value 
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Figure 4. We show, for a typical event, the time evolution of the specific heat per 
particle at constant volume Cy for N = 500 and U = 0.6. Cy assumes negative values 
in the non-equilibrium regime (up to time t ~ 3000). The dashed straight line gives the 
equilibrium value. 



smaller than the canonical prediction. The same behavior appears in a whole energy 



region 0.5 < U < Uc and, as shown in ref. |]T9|,^, it corresponds to the presence of Quasi- 
Stationary States (QSS). The plateaus are longer the greater the system size, and in the 
A^ — >■ oo the system does not relax to the standard canonical equilibrium, remaining 
forever in the QSS [^. The QSS represent a different dynamical equilibrium for non- 
extensive systems, and according to an intriguing scenario recently proposed by Tsallis 
[EHl, the reason could be the inversion of the time limit with the size limit, 



an inversion 



which is implicit in our numerical simulations. 

In fig. 3 we report the same quantities shown in fig.l, now calculated for a very short 
time, i.e. at t = 3000 instead of t = 2 ■ 10^. For this time, average variables, such as T 
and M, have almost reached the canonical equilibrium for a large range of initial energies 
U, see panel (a) and (b). There are, however, relevant differences with respect to the 
equilibrium values in the region 0.5 < U < Uc, in correspondence to the presence of the 
QSS, where the caloric curve T(U) shows a backbending typical of a first-order phase 
transition, see in particular the magnification reported in the inset of panel (a). On the 
other hand, fluctuations are affected by non-equilibrium effects in the whole spectrum of 
energy. In fig. 3(c) the kinetic energy fluctuations are bigger than the equilibrium ones. 
As a consequence, see eq.4, the specific heat per particle reported in fig.3(d) assumes 
negative values and the behavior of Cy versus U thus obtained is similar to the one 
found in nuclear multifragmentation data, see fig. 4 of ref.[|l3l. In particular the authors 



of ref. ||13[ claim that such a negative branch in the heat capacity is a direct evidence of 
a first-order phase transition. Here we have a practical example that also a system with 
a second-order phase transition can show a negative specific heat and a backbending in 
the caloric curve, due to the non perfect equilibration. This effect increases with the size 
of the system and becomes stable in the limit N —>■ oo |2^j2^. This dynamical effect 
that simulates a first-order phase transition can be somehow explained by the presence 
of superdiffusion and Levy walks in the out-of-equilibrium regime |]19i, that implies the 



coexistence of a liquid (clustered particles) and a gas (free particles) phase ||26[| . 

It is interesting to investigate the relaxation of the specific heat to the equilibrium 
value. In fig.4 we report Cy, calculated by means of eq. §, versus time for the case 
U=0.6. Cv is negative in the out-of-equilibrium regime, in correspondence of QSS, and 
converges very slowly to the equilibrium value. It is therefore necessary to wait for very 
long times, longer than for the temperature thermalization, to have also the thermalization 
of specific heat to positive values. A more complete study of the specific heat in HMF is 
in preparation |p^6||. We notice in conclusions that the slow relaxation around the critical 
point persist notwithstanding the strong chaoticity found for the microscopic dynamics 
and reported in our previous work p7| , P^|20| . The reason for that is not completely 
understood, though it has been found that the Lyapunov exponent is proportional to the 
kinetic fluctuations and some heuristic conjectures have been proposed [^ . 



4. CONCLUSIONS 

We have presented equilibrium and out-of equilibrium microcanonical simulations of 
HMF, a simple system of interacting spins. Though the model is not explicitly built to 
describe nuclear system, we think it can give interesting information for the multifrag- 
mentation phase transition. In fact, HMF can be solved in the canonical ensemble and 
has a second-order phase transition in the thermodynamic limit. Moreover, the advantage 
with respect to Ising, lattice gas models and percolation, is that HMF allows a dynamical 
microcanonical approach to the study of phase transitions in a finite system. HMF has 
a very rich dynamical behavior in a transient out-of-equilibrium regime whose timescale 
depends on the energy and on the size of the system. When the microcanonical simula- 
tions are started in an out-of-equilibrium initial state, for example in a "water bag" , we 
find the appearance of quasi-stationary-states. In correspondence we have a caloric curve 
with a well defined backbending and negative specific heat. These indications of a first- 
order phase transition, at variance with the second-order phase transition predicted in the 
thermodynamic limit, are effects of the non-equilibrium. Certainly the non-equilibrium 
features of the model are strictly linked to the long-range nature of the interaction, how- 
ever the general validity of these results is not completely clear and further work in this 
direction is needed. 

We thank X. Campi, M. Pettini and C Tsallis for useful discussions. Part of this work 
has been done in collaboration with S. Ruffo. 

REFERENCES 



1. J. Pochodzalla et al. , Phys. Rev. Lett , 75 (1995) 1040; See also the contributions 
presented at this conference by W. Miiller and U. Lynen. 

2. Y.-G Ma et ah, Phys. Lett. 5, 390 (1997) 41 and refs. therein. 

3. A. Chbihi, O. Shapiro, S. Salou and D.H.E. Gross, Eur. Phys. J. A, 5 (1999) 251. 

4. Gilkes et al, Phys. Rev. Lett. , 73 (1994) 1590; Elhott et ah, Phys. Rev. C , 49 (1994) 
3185; 

5. R.P. Sharemberg et al. submitted to Phys. Rev. Cand contribution presented at this 
conference by N.T. Porile. 

6. M. D'Agostino et al. , Nucl. Phys. A , 650 (1999) 329. 

7. D.H.E. Gross, Phys. Rep. , 279 (1997) 119 and refs therein; see also the contribution 
presented at this conference. 

8. J. Bondorf et al., Phys. Rep. , 257 (1995) 133 and refs therein. 

9. V. Latora, M. Belkacem, A. Bonasera, Phys. Rev. Lett.., 73 (1994) 1765; M. Belkacem, 
V. Latora, A. Bonasera, Phys. Rev. C, 52 (1995) 271. 

10. A. Bonasera, V. Latora and A. Rapisarda, Phys. Rev. Lett., 75 (1995) 3434 . 

11. A. Atalmi, M. Baldo, G.F. Burgio and A. Rapisarda, Phys. Rev. C, 53(1996) 2556; 
and Phys. Rev. C, 58 (1998) 2238. 

12. X. Campi and H. Krivine, Nucl. Phys. A 620 (1997) 46; X. Canipi and H. Krivine, A. 
Puente, Physica A 262 (1998) 328. see also the contribution presented by X. Campi 
at this conference. 

13. M. D'Agostino et al., Phys. Lett. B, 473 (2000) 219; see also Ph. Chomaz and F. 
Gulminelh, Nucl. Phys. A 647 (1999) 153. 

14. W. Bauer, Phys. Rev. C , 38 (1988) 1297; W. Bauer and Botvina, Phys. Rev. C, 55 
(1997) 546; 

15. L. Moretto and G.J. Wozniak, Ann. Rev. Nucl. Part. Sci. , 43 (1993) 379; see also L. 
Moretto's contribution at this conference. 

16. S. Ruffo, Transport and Plasma Physics, S. Benkadda, Y. Elskens and F. Doveil eds.. 
World Scientific, Singapore (1994), p. 114; M. Antoni and S. Ruffo, Phys. Rev. E 52 
(1995) 2361. 

17. V. Latora, A. Rapisarda and S. Ruffo, Phys. Rev. Lett. 80 (1998) 692. 

18. V. Latora, A. Rapisarda and S. Ruffo, Physica D 131 (1999) 38. 

19. V. Latora, A. Rapisarda and S. Ruffo, Phys. Rev. Lett. 83 (1999) 2104 and Physica 
A 280 (2000) 81. 

20. V. Latora, A. Rapisarda and S. Ruffo, Progr. Theor. Phys. Suppl. 139 (2000)in press, 
and ||cond-mat/00010iq| . 

21. M. Antoni, H. Hinrichsen, and S. Ruffo, to be published in Chaos solitons and fractals, 
[|cond-mat/98 100481] . 



22. T. Dauxois, P. Holdsworth and S. Ruffo, ||cond-mat/0002005|| 



23. H. Yoshida, Phys. Lett. A 150 (1990) 262. 

24. J.L. Lebowitz, J.K. Percus and L. Verlet Phys. Rev. A 153 (1967) 250. 

25. E.M. Pearson, T. Halicioglu and W.A. Tiller Phys. Rev. A 32 (1985) 3030. 

26. V. Latora and A. Rapisarda, in preparation. 

27. V. Latora and A. Rapisarda, to be published in Chaos solitons and fractals 
mat/0006TT^ . 



3onc 



28. C. Tsallis, in " Nonextensive Statistical Mechanics and Thermodynamics", eds. S.R.A. 



Salinas and C Tsallis, Braz. J. Phys. 
Today (2000) in press. 



29, 1 (1999) ||cond-mat/9903356|| and Physics 



